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We study a simple modification of the optimized random phase approximation (ORPA) aimed at 
improving the performance of the theory for interactions with a narrow attractive well by taking into 
account contributions to the direct correlation function non-linear in the interaction. The theory is 
applied to a hard-core Yukawa and a square-well potential. Results for the equation of state, the 
correlations, and the critical point have been obtained for attractions of several ranges, and compared 
with Monte Carlo simulations. When the attractive interaction is narrow, the modified ORPA 
significantly improves over the plain one, especially with regard to the consistency between different 
routes to the thermodynamics, the two-body correlation function, and the critical temperature. 
However, while the spinodal curve of the modified theory is accessible, the liquid-vapor coexistence 
curve is not. A possible strategy to overcome this drawback is suggested. 

I. THE CLOSURE 

In the last years there has been a considerable interest in model fluids in which a hard-core or very steep repulsion 
in the interparticle potential is followed by a narrow and deep attractive well. This interest mainly stems from 
the fact that such potentials provide a modelization of the interaction in many colloidal suspensions and protein 
solutions which, albeit crude, succeeds nevertheless in capturing the most important features of thermodynamics [EJ, 
structure j|, and phase behavior ||^] of these systems. When describing the properties of fluids with narrow attractive 
potentials, it appears that the second virial coefficient plays an important role. In particular, while the critical 
temperature T c depends very sensitively on the interaction range, the second virial coefficient B2 evaluated at T c 
changes quite slowly [pj. Remarkably, this remains true even when considering different interaction profiles, e.g. 
switching from a square- well to an attractive Yukawa potential. On the basis of simulation results, B2(T C ) is found 
to lie in the neighborhood of B2(T C ) ~ — 6so U> where sq = ira 3 /6 is the volume of the particle, represented as a 
hard sphere of diameter a. Other properties of fluids with short-range attractions exhibit a considerable degree of 
universality when described in terms of £?2- For instance, the range of the potential itself can be defined as that of 
an "equivalent" square- well fluid with the same value of B2 as the original interaction Q. With this prescription, 
the disappearance of a stable liquid-vapor transition characteristic of short-ranged attractive potentials |3||l|,0-f§ has 
been found to occurr at about the same range for several different interactions. It has also been observed that B2 is 
a relevant parameter in protein crystallization, the optimal conditions for crystallization corresponding to a narrow 
range of values of B2 [[H)| . 

Probably the most widespread method to model narrow attractive interactions in colloidal systems makes use of 
Baxter's exact solution of the adhesive hard sphere (AHS) model jll| within the Percus-Yevick (PY) integral equation. 
As it is well known, this can be thought of as a square- well fluid in the limit of vanishing range and infinite well depth 
so that B2 will stay finite. Therfore, the relevance of the second virial coefficient in this model stems from the fact that 
this is the quantity which allows one to make contact with the real system. It is worthwhile recalling that the existence 
of an equilibrium solution for the AHS model is a consequence of the approximate character of the PY theory. In fact, 
it has been pointed out that the AHS model is in itself thermodynamically unstable. While Baxter's solution is 
successful in describing many properties of systems with very short-ranged attractive potentials, one would also like 
to describe these systems without resorting to the sticky limit. It is well known that the most popular perturbative 
approaches in liquid-state theory, such as the mean spherical approximation (MSA) and the optimized random phase 
approximation (ORPA) pi are not suited for short-ranged interactions. The reason is that both of these theories 
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assume that the direct correlation function c(r) outside the repulsive core is simply proportional to the interaction 
w(r). This is a sound assumption whenever w(r) is long-ranged, but becomes less and less so as the range of w(r) 
decreases to a small fraction of the size of the particles. An accurate description of the correlations is provided by 
other in teg ral-equation approaches, the most successful of which is probably the modified hypernetted chain (MHNC) 
theory J14[ , in which c(r) depends on the interaction both directly and via the distribution function g(r). On the 
other hand, determining the critical point by these more sophisticated methods proves difficult or even impossible, 
while this is accomplished quite easily within the simpler perturbative schemes. Therefore, it can be interesting to 
seek for a modification of the ORPA which mantains the simplicity of the original formulation, while accounting 
more accurately for the effects on c(r) non-linear in the interaction. This could prove useful also in the context of 
more complex theories, namely the renormalization-group based hierarchical reference theory (HRT) []l5f and the 
self-consistent Ornstein-Zernike approximation (SCOZA) pq |, both of which adopt a simple ORPA-like functional 
form for c(r). These theories deliver the best results for the critical point and the liquid- vapor coexistence curve of 
simple fluids, but the inadequacy of ORPA for very short-ranged interactions partly affects their accuracy in this 
regime, especially as far as the correlations are concerned [p~7| ] . 

A possible way to modify the ORPA is suggested by the aforementioned observation that the second virial coefficient 
evaluated at the critical temperature depends little on the specific form of the interaction. For a potential v(r) which is 
the sum of a hard-core contribution and an attractive tail w(r), this behavior could be accounted for, at the mean-field 
level, by replacing the usual random phase approximation (RPA) for the direct correlation function c(r) 

c(r) = cas(r) - (3w{r) (1) 

by the expression 

c(r)=c HS (r) + e-^ p )-l > (2) 

where (3 — I/ZcbT is the inverse temperature and chsM is the direct correlation function of the hard-sphere reference 
fluid. The latter can be considered as a known quantity, for instance by resorting to the accurate Verlet-Weis 
parameterization Jig] , Eq. (|J) corresponds to a mean-field or van der Waals-like expression for the Helmholtz free 
energy per unit volume a: 

(3a = (3a HS + \p 2 J d 3 r [1 - e"^] , (3) 

where ans is the Helmholtz free energy per unit volume of the hard-sphere reference fluid. As is readily seen from 
this expression, the mean-field critical density p c does not depend on the attractive interaction w{r) and is given by 
p c = 0.249er~ 3 if the Carnahan-Starling expression for the reduced compressibility of the hard-sphere gas 
is used. As a consequence, Eq. (||) gives a universal value for B2{T C ) at the critical point: 

B 2 (T C ) = \ fd 3 r [1 - e"W f )] = |™ 3 - - * = -6.601 s , (4) 

where j3 c is the inverse critical temperature. This is not far from the above-mentioned value B^iTc) ~ — 6so reported 
in Ref. fjj. We have then considered a similar modification of the ORPA by replacing the perturbation —f3w(r) 
with the Mayer function e~@ w ^ — 1. We will refer to this approximation as non-linear ORPA. The Ornstein-Zernike 
equation which relates the direct correlation function c(r) to the radial distribution function g(r) is then closed by 
means of the following ansatz: 

g(r)=0 r<a : 

c(r) = chsO*) +e-^M - 1 r > a . { > 

This closure differs from Eq. (|2j) insomuch as here, as in the usual ORPA, we require g(r) to vanish inside the core 
because of the singular repulsion. Instead of Eq. (g) we have then: 

c(r) = c ns (r)+e-^-l + g(r), (6) 

where the function Q{r) vanishes for r > a, and is determined so that the core condition g(r) — for r < a is satisfied. 
The latter can be written via the Ornstein-Zernike equation in terms of the direct correlation function in Fourier 
space c(k): 



d 3 k lk . r c(k) 



■ e 



(2tt) 3 1 - pc{k) 



= -1, r<a. (7) 
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As is customary in ORPA calculations, we have implemented this condition in an approximate fashion by representing 
G(r) for r < a on a finite set of basis functions as a polynomial of given (typically, fourth) degree: 



g(r)=J2u j r j - 1 , r<a. (8) 

The state-dependent coefficients Uj, j — 1 . . . n, are determined by projecting Eq. (Q) on each basis function in the 
interval < r < a. If we set: 

and we indicate by f(r) the Mayer function e~ /3w ^ r ' ) — 1, Eq. (fjj) gives the following set of n non-linear equations for 
the n coefficients uf 

Pi(k) = : J - : T = -P ; (fc = 0), 1 = 1. ..n, (10) 



(2tt) 3 



1-p 



2Hs(fe)+/(*)+ELi«i^(A) 



where the hat has been used to denote Fourier transform. This system has been solved numerically by means of a 
standard Newton-Raphson iterative algorithm. Particular care has been taken to ensure convergence of the algorithm 
even for states at which the compressibility of the system is large, i.e. such that the denominator of the integrand 
in Eqs. (0), (|Io|) becomes small near k — 0. By evaluating analytically the contribution to the integral at small 
wavevector fc, we could easily achieve convergence for reduced compressibilities as large as ~ 10 8 . This enabled us 
to accurately locate the compressibility-route critical point and the spinodal curve (that is, the locus of diverging 
compressibility) predicted by the theory. In order to check the reliability of the approximate procedure adopted to 
solve Eq. (Q), we used the same method to solve the MSA for a hard-sphere fluid with an attractive Yukawa tail and 
compared the result with the exact solution of the MSA for this particular potential jl9| . The MSA is obtained from 
Eq. (^|) by using for chsM the solution of the PY equation for the hard-sphere fluid and by replacing the Mayer 
function with the perturbation —j3w(r). The comparison shows that the results obtained by the procedure described 
above are nearly undistinguishable from those of the exact solution for all the interaction ranges investigated here. 
The relative error on the reduced compressibility x ro d is typically of the order of 10~ 5 and does not exceed few percents 
for Xrcd in the range 10 4 -10 5 , while critical temperatures coincide to within four significant figures. 



II. RESULTS 

We have assessed the performance of the non-linear ORPA of Eq. (|^) by comparing its predictions for the thermo- 
dynamics, the correlations, and the critical point with those of other liquid-state theories as well as simulations. Most 
of the results were obtained for a hard-core Yukawa fluid at several values of the interaction range. The potential has 
the form 

v(r) = {_ eae - z{r -„ )/r r>(Ji (11) 

where z is the inverse-range parameter. We also performed some calculations for a hard-core square-well potential 

oo r < a , 

v (r) = { -e cr < r < (I + 5) a, (12) 
r> (l + 5)o-, 

where 5 measures the width of the attractive well. It is worthwhile observing that for this particular potential the 
closure (||) amounts to evaluating the plain ORPA at an "effective" inverse temperature /3 e ff given by /3 e ff£ = e' 36 — 1- 
This corresponds to rescaling the amplitude of the perturbation —(3w(r) so that the second virial coefficient is correctly 
reproduced. Such a prescription was in fact adopted several years ago to improve the performance of the MSA for 
the very narrow square- well potentials used to model the structure factor of microemulsions ■ We also recall that 
the above interactions can be used to mimic a Lennard- Jones (LJ) fluid by suitably adjusting their ranges. Adequate 
values for this purpose are z = 1.8cr _1 [2l| for the hard-core Yukawa potential ( |TT| ) and 5 = 0.5 for the square-well 
potential (|l2|). In the following we will deal with substantially shorter-ranged interactions. 
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In figures we report the results for the equation of state of the hard-core Yukawa fluid. Reduced temperatures 
T* = ksT/e and densities p* — pa 3 will be used throughout. We have considered both the linear and the non-linear 



ORPA, and have compared them with the MHNC |17|] and Monte Carlo (MC) simulations [|17|j22|. In order to deter- 
mine the degree of thermodynamic consistency, results were obtained by all the three usual routes to thermodynamics, 
namely the compressibility, the virial, and the energy route. We recall that these are different prescriptions to obtain 
the thermodynamics from the two-body correlations. Specifically, in the compressibility route the key to the thermo- 
dynamics is the reduced compressibility Xrcd determined as the structure factor S(k) at zero wavevector k. In terms 
of the direct correlation function c(r), Xrod is then given by 



1 f°° 

= 1 — 4-7T/9 / dr r 2 c(r) . 

Xrcd Jo 



(13) 



In the internal energy route, the starting point to determine the thermodynamics is instead the excess internal energy 
per particle u, obtained by the integral of the interaction weighted by the radial distribution function g(r): 

poo 

u = 2np / drr 2 g(r)w(r). (14) 

J a 

Finally in the virial route the thermodynamics is obtained from the pressure P as given by the virial theorem: 



0P 
P 
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a'mrr-) i I drr 3 g{r)^(r) 



(15) 



where the contact value of the radial distribution function g(a + ) appears because of the discontinuity of the hard- 
core plus tail interaction considered here. While in an exact treatment Eqs. (|l3|), (fl4|), (|l5|) are obviously bound to 
give equivalent results, this is in general not true as soon as approximations are introduced. The extent to which 
the equations of state determined via the various routes differ from one another is a useful test of the reliability of 
the theory. Figures show that the non-linear ORPA is much more consistent than the plain ORPA, especially 
for short-ranged interactions. It is known |Q that for the plain ORPA the internal energy route delivers the most 
reliable results for the equation of state. This is true for the non-linear ORPA as well. The internal energy equations 
of state are actually very similar in the two cases, although for short-ranged potentials and at low temperature 
(see figure ||) the non-linear theory appears to be slightly superior. It is worthwhile observing that, while in the 
linear ORPA the pressure P cn obtained via the internal energy route always lies between the compressibility pressure 
-Pcomp and the virial one P vlr with P V j r < P cn < F comp , this is not always the case with the non-linear ORPA. For 
instance, for an inverse interaction range z = 9cr _1 and reduced temperature T* = 0.8 the non-linear ORPA has 
again P vir < P cn < P com p (see figure ||), while for the same value of z and T* = 0.45 it gives ? comp < P cn < P V1I (see 
figure H). For short-ranged interactions the energy- route equation of state is more accurate than the MHNC, which 
underestimates the simulation data. 

The radial distribution function g(r) and the direct correlation function c(r) of the hard-core Yukawa fluid for several 
thermodynamic states and interaction ranges are shown in figures Again, the results obtained by the linear and 
non-linear ORPA have been compared with the predictions of MHNC, which is known to describe two-body correlations 
quite accurately. For g(r), MC simulation data [|l7j have also been reported. As anticipated, the plain ORPA appears 
to be inadequate to describe the correlations for short-ranged interactions. In particular, the contact value of g(r) is 
severely underestimated, with a relative error easily reaching 100%, and overall the profile of g(r) does not show enough 
structure. The non-linear ORPA performs considerably better by comparison, although some discrepancies with the 
MC data emerge, especially at low temperature and high density. These go in the opposite direction with respect to 
the linear ORPA, the contact value being overestimated, and the structure being somewhat overemphasized. As a 
consequence, as shown in the figures, the exact g(r) lies between the linear and non-linear ORPA results. Hence one 
expects that an accurate representation of the correlations can be obtained by performing a weighted average of the 
two closures. Such an approach has in fact been proposed by other authors fl24f and referred to as tail interpolation 
method (TIM). For a hard-core plus tail potential, this consists in writing the contribution of the tail interaction 
to the direct correlation function outside the repulsive core as c(r) - chsM = a[e- pw ^> - 1] + (1 - a)[~(3w(r)}, 
where the parameter a is determined as a function of temperature and density by imposing consistency between the 
virial and the compressibility routes. In Ref. |24|] this method was applied to the hard-core Yukawa potential with 
z = 1.8(T _1 as well as the Lennard- Jones potential via the Weeks-Chandler-Andersen ^5| splitting of the interaction, 
and very accurate results for g(r) were obtained. It could then be interesting to apply TIM also to narrower potentials. 
The possibility of improving the performance of the closure studied here by embedding it into a thermodynamically 
self-consistent scheme is discussed at the end of this section. 
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We observe that, unlike the linear ORPA, the non- linear closure (||) becomes exact at low density p, as it yields the 
correct zeroth-order term of the expansion of c(r) in powers of p, thereby giving c(r) ~ e -0 v ( r ) _ i 5 g( r ) e -Pv(r) f or 
p — > 0. This also implies that the second virial coefficient is correctly reproduced. Actually, it has been observed |26| ] 
that for short-ranged interactions the form g(r) ~ e~^ v ^ remains reasonably accurate near contact in a remarkably 
wide density interval. This is true in particular in the regime in which the fluid-solid equilibrium curve broadens 
up and pre-empts the gas-liquid transition. Figure Q shows the contact values of g(r) as predicted by non-linear 
ORPA, MHNC, and MC simulations for the fluid with z — 9a~ x and T* = 0.4 as a function of density. According 
to simulation studies |5), the states shown in the figure are well inside the fluid-solid coexistence region. The contact 
value of g(r) given by the non-linear ORPA is indeed very close to the value e _1 / T ~ 12.18 for densities as high as 
p* = 0.8. The insensitivity of the peak of g(r) to the density is due to the fact that the clustering of the particles 
near contact is triggered by the nearly sticky attractive interaction [ p6| . This sharply contrasts the situation found in 
Lennard- Jones- like systems away from the critical point, where the structure is mainly determined by the excluded 
volume effect, hence being essentially that of a hard-sphere fluid at the same density. 

We now turn our attention to the critical point. We have obtained critical temperatures and densities both for the 
hard-core Yukawa potential of Eq. ([ll]) and the square- well potential of Eq. (^2|). The results obtained for several 
interaction ranges by the compressibility and the internal energy route of the linear and non-linear ORPA have been 
compared with simulation results fip^|,p^,pl|j7p|Jl^| in tables | and || 

The longest ranges considered correspond 
to the above-mentioned LJ-like values z = 1.8<r -1 , S — 0.5 for the Yukawa and the square-well fluid respectively. 
The higher degree of thermodynamic consistency of the non-linear ORPA is evident from the better agreement 
between energy- and compressibility-route critical temperatures with respect to the linear ORPA. In particular, the 
energy-route critical temperature of the non-linear ORPA is in quite good agreement with the simulation results 
in the whole range interval considered, although a precise assessment of the accuracy of the theory in this respect 
is prevented by the rather large uncertainties that usually affect critical point estimates based on simulations, as 
is apparent from the discrepancies between the simulation results reported in table |n|. The most reliable results 
are obtained by supplementing simulations with a finite-size scaling analysis; however, to our knowledge this was 
performed only for the longest ranges among those reported in the tables [ |16|j2ljj3l|| . It is clear in any case that the 
compressibility route of linear ORPA gives the worst predictions, by severely underestimating the critical temperature 
as the interaction range decreases. For the critical density p c the non-linear ORPA performs less accurately than for 
T c : the compressibility route overestimates p c with respect to simulation, while the energy route underestimates it. 
According to what observed above, for the square-well fluid the compressibility-route reduced critical temperatures 
T*j and T* nl predicted respectively by the linear and non-linear ORPA are linked by the relation 1/T*, = e 1 ^^ — 1. 
However, this relation does not hold for the critical temperatures given by the internal energy route. 

Since in introducing the closure (|^) we were inspired by the empirically observed behavior of the second virial 
coefficient at the critical temperature, it is worthwhile considering the results for B2(T C ) given by the theory. These 
have also been reported in the tables. Moreover, in figure ||we have plotted the reciprocal of the critical temperature of 
the square- well fluid as a function of the logarithm of the inverse well width 1/5. If B 2 (T C ) were rigorously independent 
of (5, this plot would yield a straight line in the small-5 limit. The slow variation of B2(T C ) with S gives a nearly 
linear dependence of 1/T C on ln(l/<5) in a relatively wide interval of 5, with the non-linear ORPA closely following 
the simulation results. However, neither the simulation data nor the theory show evidence of B2(T C ) saturating to a 
constant value in the interval considered. We observe that the linear ORPA gives rather poor results for B2(T C ). In 
particular, the values obtained via the compressibility route of linear ORPA are way too large in modulus and even 
show the wrong trend because of the strong underestimation of T c . 

Below the critical temperature of the compressibility route, both the linear and the non-linear ORPA give a 
spinodal curve, along which the compressibility given by Eq. (|r|) diverges. Inside the region bounded by this curve 
the approximation cannot be solved, as the integrand in the l.h.s. of Eq. ( |lC)| ) would develope a pole at non-zero 
wavevector and the integral would diverge. In order to obtain the coexistence curve, the problem arises of avoiding 
this forbidden region when determining the pressure and the chemical potential needed to impose the thermodynamic 
equilibrium conditions. In the linear ORPA, this is achieved via the internal energy route, namely by getting the excess 
Hclmholtz free energy by integrating along each isochore the excess internal energy given by Eq. (flj). The success 
of this procedure depends on the fact that, as shown in tables |, ||, the critical temperature of the compressibility 
route of linear ORPA is considerably lower than that of the internal energy route, at least for the simple model fluids 
considered here. Therefore, the forbidden region lies well inside the internal energy coexistence boundary. However, 
we have found that this method is not effective when dealing with the non-linear ORPA. In fact, as observed above, 
in this case the difference between the critical temperatures given by the compressibility and the internal energy route 
is relatively small. As a consequence, the relative location of the energy-route critical point and the spinodal curve is 
such that a substantial portion of the coexistence curve should lie inside the forbidden region. This prevented us from 
determining the coexistence curve, save for a tiny part corresponding to the very top. The situation just described 
is illustrated in figure 0, which shows the spinodal curve predicted by the non-linear ORPA for a hard-core Yukawa 
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fluid with z = 4er~ 1 , together with the locus of points at which the internal-energy route compressibility diverges. For 
densities just above the critical one, the latter curve bumps into the spinodal. This rapidly widens as the temperature 
is lowered, so that also the high-density branch of the coexistence curve is swallowed by the forbidden region. Such a 
feature is clearly quite a serious limitation of the non-linear ORPA, especially since the thermodynamic consistency of 
the theory is not so high as to allow one to reliably determine the pressure and the chemical potential along a mixed 
integration path in the temperature-density plane, as it is often done with other integral-equation theories |p^ , ^o| . 
A possible strategy to overcome it consists in using a functional form of c(r) similar to that considered here in the 
context of a more sophisticated theory, such as the SCOZA or the HRT mentioned in the previous section. This could 
be accomplished by multiplying the Mayer function in Eqs. (||), (||) by a state-dependent amplitude A(p, T). In the 
SCOZA A(p, T) is determined by imposing consistency between the compressibility and the energy route, while in 
the HRT one requires the compressibility rule ( p"3| ) to be satisfied along a process in which the Helmholtz free energy 
of the system is obtained by gradually switching on the long-wavelength Fourier components of the interaction. Both 
of these procedures are expected to lead to a well-defined coexistence curve. In fact, in the SCOZA the spinodal 
curve and the curve along which the internal-energy route compressibility diverges are forced to coincide. Since the 
coexistence boundary lies outside the latter curve, it will lie outside the forbidden region as well. On the other 
hand the HRT, thanks to its renormalizaton-group structure, preserves the convexity of the free energy for every 
thermodynamic state. No forbidden region appears, and the two-phase domain is correctly obtained as the locus of 
flat pressure-density isotherms p5[ . As stated above, the accuracy of both the SCOZA and the HRT for short-ranged 
interactions is affected by the functional form adopted for c(r), similar to that of the linear ORPA. The predictions 
for c(r) given by the linear and non-linear ORPA for a number of different states of the Yukawa fluid have already 
been compared to the accurate MHNC results in figures [|-^| In figure [lO] the direct correlation function outside the 
repulsive core for the Yukawa fluid with z = 9<r _1 , T* = 0.4, p* = 0.8 shown in figure^ has been rescaled by its value 
at contact c(er + ). It appears that, even when both the linear and the non- linear ORPA differ considerably from the 
MHNC as for the state considered here, the overall profile of c(r) given by the non-linear ORPA can be adjusted so 
as to reproduce that of the MHNC fairly well by suitably rescaling its amplitude, while this is not the case with the 
linear theory. Accordingly, we expect that implementing a closure similar to that of Eq. (^|) within the HRT or the 
SCOZA scheme along the lines sketched above will lead to an accurate description of thermodynamics, correlations, 
and phase diagram of fluids with narrow attractive potentials. We intend to pursue this investigation in the future. 

III. CONCLUSIONS 

We have considered a simple modification of the ORPA aimed at improving the performance of this theory for 
interactions which present a narrow and deep attractive well. Its predictions have been compared with those of linear 
ORPA and MHNC as well as simulation data. The equation of state is accurately reproduced, and the degree of 
thermodynamic consistency is considerably higher than in the plain ORPA. Also the two-body correlations are more 
accurate than those of plain ORPA, although overall less accurate than in MHNC. The critical temperatures agree 
well with simulation results. In particular, the second virial coefficient at the critical temperature B2(T C ) is found to 
depend rather weakly on the interaction range. Below the critical temperature the spinodal curve predicted by the 
theory can be easily obtained, but unfortunately most of the liquid- vapor coexistence curve is not accessible, as it lies 
inside the region where no solutions can be found. This drawback could be overcome by using a functional form of 
the direct correlation function similar to that considered here in the context of a more sophisticated approach such 
as the HRT or the SCOZA. 
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-6.14 


0.576 
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-5.92 
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0.412 


0.50 


-5.51 


ORPA (comp) 


1.037 


0.312 


-7.89 


0.542 


0.360 


-11.10 


0.387 


0.392 


-15.35 


0.213 


0.460 


-46.70 


ORPA (energy) 


1.243 


0.318 


-5.61 


0.771 


0.374 


-5.16 


0.624 


0.420 


-4.75 


0.451 


0.536 


-4.05 


ORPANL (comp) 


1.221 


0.326 


-5.82 


0.749 


0.388 


-5.52 


0.602 


0.432 


-5.29 


0.430 


0.504 


-4.77 


ORPANL (energy) 


1.210 


0.270 


-5.92 


0.736 


0.302 


-5.75 


0.591 


0.324 


-5.55 


0.424 


0.362 


-4.98 



TABLE I. Critical constants of the hard-core Yukawa fluid. The reduced critical density p* and temperature T* from the 
compressibility and energy route of linear and non-linear ORPA (indicated as ORPANL) are compared with simulation results. 
The second virial coefficient at the critical temperature B 2 (T C ) is in units of the hard-sphere volume 7r<7 3 /6. MCI: Monte Carlo 
simulation data with finite-size scaling from Ref. Jl6]] . MC2: Monte Carlo simulation data from Ref j?]]. MC3: Monte carlo 
simulation data from Ref. M. 
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-5.93 
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0.881 


0.348 


-16.05 


0.624 


0.388 


-21.40 


0.423 


0.420 


-32.73 


0.158 


0.537 


-737.88 


ORPA (energy) 


1.327 


0.304 


-6.68 


1.083 


0.354 


-5.72 


0.846 


0.438 


-4.61 


0.529 


0.686 


-3.43 


ORPANL (comp) 


1.319 


0.348 


-6.78 


1.045 


0.388 


-6.26 


0.824 


0.420 


-5.01 


0.502 


0.537 


-4.38 


ORPANL (energy) 


1.282 


0.228 


-7.22 


1.038 


0.266 


-6.36 


0.806 


0.330 


-5.37 


0.497 


0.450 


-4.57 



TABLE II. Same as table I for the square-well fluid. MCI: Monte Carlo simulation data with finite-size scaling from Ref. J31| . 
MC2: Monte Carlo simulation data from Ref. | p7[ . MC3: Monte Carlo simulation data from Ref. [Q. MD: molecular dynamics 
simulation data from Ref. |28|. 
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FIG. 1. Compressibility factor of the hard-core Yukawa fluid (see Eq. (]ll|)) with inverse range z = 4<r~ at reduced tem- 
perature T* = 0.7, corresponding to T/T c ~ 1.2. The compressibility, virial, and energy routes of the linear ORPA (upper 
panel) and non-linear ORPA (lower panel) are compared with MHNC and MC results. Upper panel: dotted line, linear ORPA 
compressibility route; dashed line, linear ORPA virial route; solid line: linear ORPA energy route; crosses, MHNC; and circles, 
MC simulation results [E2|. Lower panel: same notation as in the upper panel for the non- linear ORPA. 



9 




10 




11 




FIG. 4. Radial distribution function g(r) (upper panel) and direct correlation function c(r) (lower panel) of the hard-core 
Yukawa fluid with z — 4<r -1 at reduced density p* = 0.7 and temperature T* = 0.62, slightly above its critical value. Solid 
line, linear ORPA; dotted line, non-linear ORPA; long-dashed line, MHNC; and crosses, MC data from Ref. Wm. 
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FIG. 5. Same as figure ^for z = 9a \ p* = 0.5, T* = 0.5. According to simulation studies this state point lies somewhat 
above the fluid-solid line. 
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FIG. 7. Upper panel: contact value g{o + ) of the radial distribution function of the hard-core Yukawa fluid with z — 1.8a 1 , 
T* = 1.5 and z = 9a" 1 , T* — 0.4. Both isotherms lie slightly above the critical temperature for the corresponding z. Lower 
panel: contact value g(u + ) rescaled by its low-density limit e~ /3v ^ a '. Triangles, non-linear ORPA; crosses: MHNC; and circles, 
MC results Jl^. The shorter-ranged interaction has a much larger contact value, which remains nearly constant from p* — 0.4 
to p* = 0.8 In the same interval, the contact value of the longer-ranged interaction increases by nearly a factor 2 on increasing 
P*- 
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FIG. 8. Inverse reduced critical temperature 1/T* of the square- well fluid as a function of the width 8 of the attractive well 
(see Eq. (|l2[)). The points shown correspond to 8 — 0.5, 8 = 0.375, 8 — 0.25, and 8 — 0.1. Crosses: linear ORPA compressibility 
route; squares, linear ORPA energy route; triangles, non-linear ORPA compressibility route; diamonds, non-linear ORPA energy 
route; circles, MC simulation data from Ref. Q; and inverted triangles, MC simulation data from Ref. pT] . Lines are a guide 
for the eye. The nearly linear dependence of 1/T* on the logarithm of the inverse range 1/8 stems from B2(T C ) being a slowly 
varying function of 8. 
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FIG. 9. Dot-dashed line, spinodal curve of the hard-core Yukawa fluid with inverse range z — 4a~ according to the 
compressibility route of non-linear ORPA. No solutions are found in the region bounded by this curve. Long-dashed line, 
locus of diverging compressibility according to the energy route of non-linear ORPA. Note that this curve bumps into the 
spinodal. Dashed line, locus of diverging compressibility according to the energy route of linear ORPA; solid line, coexistence 
curve according to the energy route of linear ORPA. The compressibility-route spinodal of linear ORPA lies well inside the 
coexistence curve and is located below the plane of the figure. 
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